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ABSTRACT 

We present new observations of the fundamental ro-vibrational CO spectrum 
of V1647 Ori, the young star whose recent outburst illuminated McNeil's Neb- 
ula. Previous spectra, acquired during outburst in 2004 February and July, had 
shown the CO emission lines to be broad and centrally peaked — similar to the 
CO spectrum of a typical classical T Tauri star. In this paper, we present CO 
spectra acquired shortly after the luminosity of the source returned to its pre- 
outburst level (2006 February) and roughly one year later (2006 December and 
2007 February). The spectrum taken in 2006 February revealed blue-shifted CO 
absorption lines superimposed on the previously observed CO emission lines. The 
projected velocity, column density, and temperature of this outflowing gas was 
30 km s-\ X 10^^ cm ^, and TOO^I^^qq K, respectively. The absorption lines 
were not observed in the 2006 December and 2007 February data, and so their 
strengths must have decreased in the interim by a factor of 9 or more. We discuss 
three mechanisms that could give rise to this unusual outflow. 
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Subject headings: accretion, accretion disks — stars: flare — stars: formation — 
stars: individual (V1647 Ori) — stars:winds, outflows 

1. Introduction 

McNeil's Nebula was recently illuminated by the outburst of V1647 Ori (McNeil 2004), 
a young star that is embedded in the Lynds 1630 dark cloud and coincides with the 850//m 
continuum source OriBsmmSS (Mitchell et al. 2001). V1647 Ori has a flat SED in the mid- 
infrared, and is thus a Class I young stellar object (YSO; Andrews et al. 2004). V1647 Ori 
underwent a similar outburst as recently as 1966 (Aspin et al. 2006), indicating that V1647 
Ori is also an EXor. Such prc-main sequence stars undergo eruptive events that dramatically 
increase their luminosity for periods of months to years (Hartmann 1998). The outbursts arc 
thought to be triggered by a rapid increase in the stellar accretion rate (Hartmann & Kenyon 
1996). The eruption in November 2003 of V1647 Ori lasted two years. During the outburst, 
the star brightened by a factor of 50 in X-rays (Kastner et al. 2006), a factor of 250 in the 
red (6 mag in the Rc band; Fedele et al. 2007a; Briceno et al. 2004), a factor of 15 in the 
near-IR(~3 mag in the J, if, and K bands; Reipurth Sz Aspin 2004), and a factor of ~15 
at wavelengths from 3.6;um to 70 /im (Muzerolle et al. 2005). From the overall brightening 
of the source, Muzerolle et al. (2005) concluded that the bolometric luminosity increased 
by a factor of 15 to 44Lq (sec also Andrews ct al. 2004) and that the stellar accretion rate 
increased from ~ 1O~^M0 yr^^ to ~ lO^^M© yr^^. Similarly, Gibb et al. (2006) inferred a 
stellar accretion rate of 3 — 6 x 1O~^M0 yr~^ from the luminosity of the Br7 emission one 
year later. This is somewhat larger than the typical accretion rate of a young low mass star 
(10~^ — IO^^Mq yr~^; Bouvier et al. 2007), yet lower than is expected for a star of the FUor 
type (~ 10-"^ Mq yr-i; Hartman & Kenyon 1996). 

During the onset of the outburst of VI 647 Ori, observations of atomic hues with P Cygni 
profiles provided evidence for a hot {T ~10,000 K), high velocity {v — —400 km s~^) wind 
(Briceno et al. 2004; Reipurth et al. 2004; Vacca et al. 2004; Walter et al. 2004; Ojha et 
al. 2006; Fedele et al. 2007a) with a mass-loss rate of Mwind=4xl0~^ M© yr~^ (Vacca et al. 
2004). This mass loss rate is much lower than that of the typical FUor (Hartmann & Kenyon 
1996) and comparable to that of a classical T Tauri star (cTTS; Hartigan ct al. 1995). The 
absorption component of the P Cygni profile of several lines (e.g. Pa/3) disappeared within 
a few months following the peak of the outburst in early 2004 (Gibb ct al. 2006). However, 
P Cygni structure in the Ha profile indicated that a weaker wind continued throughout the 
outburst phase (Ojha et al. 2006; Fedele et al. 2007a). 

In contrast to the hydrogen and helium lines, the fundamental near-infrared ro- vibrational 
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emission lines of CO, observed 2004 February 27, were broad, centrally peaked, and compat- 
ible in their intensity with an excitation temperature of 2500 K (Rettig et al. 2005). The 
width of the lines was shown to be consistent with Keplerian orbital motion of the gas within 
the inner disk surrounding the central star, similar to the broad emission line profiles that 
are observed around cTTSs and Herbig Ac/Be stars (HAeBes; Najita et al. 2003; Blake & 
Boogert 2004). A later observation, obtained on 2004 July 30, showed that the CO lines 
remained broad but the temperature of the gas decreased to 1700 K (Gibb et al. 2006). 
Neither observation showed any indication of CO in an outflow, as a blue shifted absorption 
component was not detected. 

We report foUowup observations of CO from V1647 Ori, which were acquired 2006 
February, 2006 December, and 2007 February. By the time of our initial 2006 observation, 
V1647 Ori had returned to quiescence and the absorption component in the Ha line profile 
had disappeared (Fedele et al. 2007a and references therein). Presumably the accretion rate 
had fallen by two orders of magnitude to its pre-outburst accretion rate as the star faded to 
its pre-outburst brightness (sec MuzeroUe et al. 2005). As suggested by the work of Najita 
et al. (2003), only a continued decrease in CO line intensity from the warm gas in the disk 
was therefore expected. However, the first post-outburst observation revealed the striking 
metamorphosis of these lines from centrally peaked emission features to emission lines with 
blue-shifted absorption. Subsequently, by late 2006 and early 2007, the CO emission fines 
returned to their original centrally peaked structure, indicating that the production of the 
outflow diminished within one year of the end of the outburst. In this paper we discuss three 
scenarios that can give rise to such a phenomenon. 



2. Observations & Data Reduction 

V1647 Ori was observed with NIRSPEC at the W. M. Keck Observatory on 2006 Febru- 
ary 17 using one M-band setting (Table 1; Figs. 1 and 2). NIRSPEC is a high-resolution 
{X/SX ~ 25, 000), near-infrared (1-5 //m), cross-dispersed spectrometer (McLean et al. 1998). 
The NIRSPEC data were reduced and cafibrated using standard reduction tecfiniques, which 
arc described in DiSanti et al. (2001), Brittain et al. (2003) and Brittain (2004). Additional 
observations of VI 647 Ori were taken witfi tfie Pfioenix spectrometer on Gemini Soutfi on 
2006 December 01, 2007 February 01, and 2007 February 06, using one setting centered at 
2011 cm~^ (Table 1). Pfioenix is a fiigfi-resolution {X/6X ~ 50,000), near-infrared (l-5;um) 
spectrometer covering 1550 km s^^ per observation (Hinkle et al. 1998; Hinkle et al. 2000; 
Hinkle et al. 2003). Tfie same reduction algoritfims used on tfie NIRSPEC data were appfied 
to tfie Pfioenix data. 



-4- 



Spectra of the P30 and P31 CO lines from the four epochs spanning 2004-2007 are 
presented in Figure 1 in order to illustrate how the CO profiles have evolved over time. 
Gaps in the spectra represent frequencies with less than 50% atmospheric transmittance. 
The sample error bars presented for each spectrum represent the standard deviation of the 
continuum. Because the individual Phoenix spectra have lower S/N, we plot just the average 
of those datasets, weighting each observation by the standard deviation of the continuum. 
The equivalent widths of the emission and absorption lines measured from the new datasets 
are presented in Table 2. Upper limits on the absorption lines in the Phoenix spectrum were 
calculated by adopting the width of the absorption lines observed in 2006 February, setting 
the depth of the line equal to the standard deviation of the continuum, and assuming a 
gaussian line profile. 



3. Results 

The barycentric velocity of L1630 is +26 km s"^ (Lada, Bally, & Stark 1991; Gibb 1999; 
Mitchell et al. 2001). Since the velocity of young embedded stars is typically within a few 
km s~^ of the surrounding cloud, we adopt this as the barycentric velocity of V1647 Ori. 
The emission lines, corrected for the motion of the earth, are centered on the barycentric 
velocity of the star under the assumption that it has the same velocity as L1630 (Fig. 1). 

Rettig et al. (2005) noted that the unblended CO emission lines observed during the 
early phase of the outburst of V1647 Ori were symmetric and centrally peaked, hence they 
concluded that the CO emission originated in the inner disk, as in cTTSs and HAeBes. 
The CO emission lines from these sources are generally centrally peaked, although a few 
of them exhibit double-peaked structure due to rotational broadening. The broadening of 
the emission lines, the absence of blue-shifted absorption, and the excitation temperature 
indicate that the emitting gas around cTTSs and HAeBes lies within a circumstellar disk 
(Najita et al. 2000). Our 2006 February observation of V1647 Ori, however, reveals the 
surprising appearance of a blue-shifted absorption feature superimposed on the CO emission 
lines. This absorption feature was blue-shifted from line center by 30 km s~^ and had a width 
of 20 km s~^. Roughly one year later, in 2006 December and 2007 February, the strength of 
the absorption decreased by at least a factor of 9 in such prominent lines as P30 (Table 2) 
and was no longer observed (Fig. 1). 

An excitation plot for the CO emission lines observed in 2006 February (after the star 
had returned to its preoutburst brightness) implies that the effective rotational temperature 
of the emitting gas was 1400±200 K. To estimate the temperature of the absorption lines, 
we fit the observations with a synthetic spectrum, which we calculate by assuming the 
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lines have a gaussian profile, the ^^CO/^^CO ratio is 60, and the gas is characterized by a 
single temperature. The free parameters in our model are the intrinsic line width, b, the 
CO column density, N(CO), and the rotational temperature of the gas, T. We present the 
predicted equivalent widths of ^^CO absorption lines in Table 2 and three models in Figure 
2. The parameters for each model are as follows. Red spectrum: T = lOOOK, b = 2.5 km 
s-\ N(CO)= 2x10^^ cm-2; green spectrum: T = TOOK, b = 2.0 km s-\ N(C0)=3 x 10^^ 
cm-2; and blue spectrum, T = 600K, b = 1.8 km s^^ N(CO)=5xlO^^ cm^^. The three sets 
of parameters do an equally good job of fitting the low-J ^^CO lines. If the gas is lOOOK 
or hotter, the high-J hues are too deep. If the gas is 600K or cooler, the low-J ^^CO hnes 
from 2127-2151 cm~^ (R9-R15) become too strong. Our la upper limit on the observed 
equivalent widths of ^'^CO lines not blended with other features is 0.006 cm~^. For gas at 
600K, TOOK and lOOOK, the predicted strengths of these lines are ~0.018 cm^^, ~0.011 cm~^ 
and ~0.006 cm^^, respectively. The model with T = TOOK (green) provides the best overall 
fit, so we conclude that T = 700tfoo K, N(CO)= 3^1 x 10^^ cm'^, and b = 2.0^0.1 km s-\ 



4. Discussion 

During quiescence, V1647 Ori is a class I YSO (Andrews et al. 2004). However, its mass- 
loss rate during outburst was similar to a strongly accreting cTTS (Vacca et al. 2004). While 
outflows from cTTSs and HAeBes are common, fundamental ro-vibrational CO emission 
lines with a blue-shifted absorption component have not been observed around any of the 
more than 300 such sources observed to date (Najita et al. 2000, 2003; Blake & Boogert 
2004; Rettig et al. 2006; Brittain et al. 2007; J. Brown private communication). Further, 
Class I YSOs such as GSS 30 IRS 1, HL Tau and RNO 91 do not show ro-vibrational CO 
emission hnes with blue-shifted absorption components (Pontoppidan et al. 2002, Brittain 
et al. 2005. and Rettig et al. 2006, respectively). In contrast to these systems, the FUor 
V1057 Cyg has shown blue-shifted overtone ro-vibrational CO absorption lines (Hartmann 
et al. 2004). While the similarity of the unusual CO outflows is intriguing, there are 
important differences between FUors such as V1057 Cyg and EXors such as V1647 Ori. 
First, FUors have greater accretion rates (~10~'* Mq yr^^) and more extreme winds than 
EXors (Hartmann & Kenyon 1996). Secondly, CO is always and on/y detected in absorption 
toward FUors, and is thought to originate in the accretion disk (Calvet et al. 1993; Calvet 
et al. 1991). Despite the significant differences between V1057 Cyg and V1647 Ori, both 
stars have undergone major eruptions that have generated outflows in CO, suggesting that 
a different mass-loss mechanism may be at work in these systems than in other young stars. 

Given the unusual nature of the outflowing CO from V1647 Ori, it is of interest to 
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consider the mechanisms that could give rise to this phenomenon. One possibihty is that 
the absorbing CO condensed out of the hot outflowing wind that was observed during the 
outburst. There were two outflow components noted in the Ha absorption feature: a variable 
component at 400 km s~^ and a steady component at 150 km s~^ (Fedele et al. 2007a). If the 
lower- velocity gas decelerated at a constant rate to 30 km s~^ (the velocity of the outflowing 
CO), by the time of our second CO observation in 2004 July the wind would have expanded 
to 10 AU. There was no evidence of blue-shifted CO absorption on this date (Fig. 1). By 
2006 February, when the absorption was observed, the moderate velocity wind would have 
expanded to nearly 40 AU. However, it seems unlikely that the wind could still retain a 
kinetic temperature of TOOK at a distance of 40AU from the star. Furthermore, it is not 
clear why CO would condense out of this outflow to reveal warm, blue shifted absorption 
but not out of any of the other outflows with similar or even greater mass-loss rates. Thus 
we conclude that this scenario is unlikely. 

A second possibility is that the CO absorption formed in a shell of material that was 
swept up by the atomic wind. In this case the absorption did not appear until the column 
density of material was sufficient to produce measurable absorption, and the heating was 
the result of the interaction of the wind with the shell. When the mass loss decreased at 
the end of the outburst phase, this heating was eventually shut down. If the CO/H2 ratio 
in such a shell was similar to that of a dense molecular cloud, 1.5 xlO~^, then the column 
density of gas was 2x10^^ cm~^. Adopting a normal interstellar extinction-to-gas ratio, 
i.e.. Ay — 5.6 X IO^^A^h mag cm^ atom~^ (Bohhn et al. 1978), we find that the observed 
column density of CO corresponded to ~20 mags of visible extinction. This is a lower hmit, 
as it is possible that selective dissociation of gas in the nebula could drive down the relative 
abundance of CO. However, the extinction measured on the line of sight toward V1647 Ori 
appears to have remained relatively unchanged over the entire course of the outburst at ~11 
mags (Vacca et al. 2004; Gibb ct al. 2006), of which 6.5 mags is due to the nebula (Fedele 
et al. 2007b). There is no evidence for an additional 10-20 mags of extinction toward V1647 
Ori, and so we conclude that this scenario is also unlikely. 

A final scenario we consider is that the CO outflow was launched in response to the 
reorganization of the stellar magnetic fleld following the sharp drop in the accretion rate. 
Pre-main sequence stars tend to have kilogauss magnetic fields which mediate stellar accre- 
tion and outfiows (e.g. Johns-KruU et al. 2000). This stellar magnetic field truncates the 
accretion disk where the ram pressure from accretion balances the magnetic pressure from 
the magnetosphere (Camenzind 1990). Consequently, the truncation radius of the disk, Rt, 
is inversely and nonlinearly proportional to the accretion rate, M, and given by Rt oc M~^/'^ 
(e.g. Bouvier et al. 2007). Thus the two order of magnitude change in the stellar accretion 
rate experienced by V1647 Ori would have resulted in the truncation radius being shifted by 
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nearly a factor of four. 

Two years into the outburst, VI 647 Ori rapidly faded by a factor of 40 in the R^-band in 
just 180 days to return to its pre-outburst level (Fedele et al. 2007a). The mid- infrared flux 
also returned to its pre-outburst level, indicating that the drop in the optical/near-infrared 
hghtcurve was due to the intrinsic fading of the source (Fedele et al. 2007a). This sharp 
drop in the luminosity of the source indicates that the accretion rate fell rapidly as the star 
returned to quiescence. In response to the drop in ram-pressure from the accretion flow, the 
truncation radius was pushed back by the magnetosphere. Evidence from simulations sug- 
gests that disk-magnetosphere systems tend to form outflows when the system is undergoing 
the greatest amount of dynamical rearrangement (e.g.. Fig. 7 in Balsara 2004). It is possible 
that the realignment of the magnetic field as it pushed out against the circumstellar disk 
resulted in a warm, shortlived outflow, an outflow that is not observed toward any other 
cTTSs or HAeBes. 

While virtually all accreting low-mass stars drive outflows, the CO outflow from V1647 
Ori is highly unusual. Indeed, the transformation of centrally peaked ro-vibrational CO 
emission lines to CO emission lines with blue-shifted absorption is unique. We suggest that 
the mechanism responsible for producing this outflow is distinct from the one that drives the 
outflow from typical cTTSs. The coincidence between the rapid fading of V1647 Ori and 
the subsequent observation of the CO outflow hints at a connection. Better sampling of the 
fundamental ro-vibrational CO spectrum of EXors as they brighten and fade is crucial for 
determining whether this coincidence is significant. The rapid and dramatic change in the 
accretion rate that characterizes the EXor phenomenon provides an important opportunity 
to study the interplay between stellar accretion, the inner disk, and outflows. This insight 
is key to reaching a satisfactory theoretical understanding of these events. 

The data presented herein were obtained (in part) at the W.M. Keck Observatory, 
which is operated as a scientiflc partnership among the Cahfornia Institute of Technology, 
the University of Cahfornia and NASA. The Observatory was made possible by the generous 
flnancial support of the W.M. Keck Foundation. Also based in part on observations obtained 
at the Gemini Observatory, which is operated by AURA, under a cooperative agreement 
with the NSF on behalf of the Gemini partnership. The Phoenix infrared spectrograph was 
developed and is operated by NOAO. The Phoenix spectra were obtained as part of program 
GS-2006A-DD-1 and GS-2006B-DD-1. S.D.B. performed this work (in part) under contract 
with JPL funded by NASA through the Michelson Fellowship Program. 
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Table 1. Log of Observations 



Date 


Telescope/ Instrument 


Spectral Grasp 


Integration 


S/N 






(cm~^) 


(seconds) 




2006 Feb 17 


KECK 2/NIRSPEC 


2151-2119, 


480 


29 






2018-1987 




28 


2006 Dec 01 


Gemini S./PHOENIX 


2016-2006 


960 


6 


2007 Feb 01 


Gemini S./PHOENIX 


2016-2006 


4800 


11 


2007 Feb 06 


Gemini S./PHOENIX 


2016-2006 


3360 


8 
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Table 2. Equivalent Widths of Observed ^^CO Lines and Model Predictions^ 



Transition Drest E.W.^ E.W.'^ E.W.'^ E.W.'^ E.W.'' 

absorption T = lOOOK T = TOOK T = 600K emission 
(cm-i) (10-2 cm-i) (10-2 cm-^) (10-^ cm-^) (10-^ cm-^) (10-^ cm-^) 



a) NIRSPEC (2006 February) 



Rl 


2150.86 


4.41 ± 0.60 


4.52 


4.69 


4.81 


(3.4) 


RO 


2147.08 


3.49 ± 0.60 


3.30 


3.98 


4.29 


(3.4) 


PI 


2139.43 


3.40 ± 1.00 


3.27 


3.96 


4.26 


(1.7) 


P2 


2135.55 


5.17 ± 0.80 


4.47 


4.64 


4.77 


(1.7) 


P3 


2131.63 


4.33 ± 0.60 


5.05 


4.97 


5.02 


(2.2) 


P4 


2127.68 


5.03 ± 0.60 


5.39 


5.18 


5.17 


-3.64±1.72 


P5 


2123.70 


6.21± 0.60 


5.62 


5.31 


5.28 


-3.59±1.72 


P6 


2119.68 


5.09± 1.00 


5.79 


5.40 


5.35 


-3.74±1.71 


P30 


2013.35 


3.57 ± 0.36 


4.45 


3.42 


3.25 


-3.80±0.39 


P31 


2008.53 


2.42 ± 0.24 


4.23 


3.15 


2.97 


-4.05±0.35 


P32 


2003.67 


1.44 ± 0.94 


4.00 


2.85 


2.66 


-3.99±0.34 


P33 


1998.78 


2.86 ± 0.64 


3.75 


2.53 


2.32 


-2.22±0.41 


b) PHOENIX (2006 December and 2007 February) 


P30 


2013.35 


(0.4) 








-2.35±0.31 


P31 


2008.53 


(0.8) 








-1.07±0.63 



^Upper limit on measured equivalent width denoted by (). 
'^Observed value. 
'^Predicted value. 
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Fig. 1. — Detailed comparison of the high-J CO lines. The v=l-0 P30 and P31 lines are 
plotted over four epochs spanning three years. An error bar has been plotted on the right 
side of each spectrum representing the standard deviation of the continuum. The spectra 
have been normalized, offset, and shifted to the rest frame of the molecular cloud in which 
V1647 Ori is embedded. 
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Fig. 2. — Comparison of the v=l-0 CO spectrum after outburst (2006 February) and syn- 
thetic models. The appearance of blue-shifted absorption lines indicates the presence of an 
outflow. We have plotted three synthetic spectra that bracket a reasonable range of param- 
eters that describe this outflowing gas. The red, green and blue spectra correspond to gas 
temperatures of lOOOK, TOOK and 600K respectively. The low temperature fits to the high-J 
lines are indistinguishable, however, the lOOOK fit (red) is too deep. The CO ice feature at 
2140 cm~^ remains unchanged which is expected in light of its interstellar origin (Rettig et 
al. 2005) and is not modeled. 



